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A  PRELIMINARY  INVESTIGATION  INTO  THE  ESTIMATION  OF 
RIVER  DEPTH  BASED  ON  MEANDER  GEOMETRY 

1.0  Introduction 

At  times  it  becomes  desirable  to  estimate  the  depths  of  rivers  using  aerial  photography  or 
other  remotely  sensed  imagery.  Due  to  water  mrbidity  in  many  of  the  cases,  optical 
techniques  using  water  color  or  laser  bathymetry  are  not  feasible.  This  preliminary  study 
investigates  the  potential  use  of  meander  geometry  based  on  empirical  equations  as  a 
possible  method  for  the  estimation  of  river  depths  for  natural  alluvial  streams  and  rivers. 

2.0  Background 

Early  studies  in  fluvial  geomorphology  undertaken  by  Luna  B.  Leopold*  and  M.G. 
Wolman  revealed  statistically  significant  correlated  relationships  between  meander  radius 
of  curvature,  wavelength,  width  and  river  depth.  These  relationships  were  found  to  be 
consistent  across  a  wide  range  of  stream  sizes;  from  small  rivulets  of  glacial  melt  to  the 
meander  structure  of  the  Gulf  Stream.  Subsequent  studies  by  Garnett  P.  Williams^  added 
additional  information  to  the  understanding  of  river  meanders  and  channel  size.  Figure  1 
illustrates  a  meander  radius  of  curvature  and  wavelength. 


WAVELENGTH  - - ». 


Figure  1.  Schematic  meander  geometry  (after  Williams  1986) 

The  following  two  figures  taken  from  several  sources  illustrate  some  of  the  relationships 
among  the  various  parameters: 


'  Leopold,  L.B.  and  Wolman,  M.G.,  1960.  River  Meander.  Bull.  Geol.  Soc.  Am.,  71:769-794 
^  Williams,  G.P.,  1986.  River  meanders  and  channel  size.  J.  Hydrol.,  88:  147-164 
Manuscript  approved  July  30, 2004. 
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Figure  2 

Relationships  with  Radius  of  Curvature  (  Rosgen,  D.,  Applied  River  Morphology  pg.  2-7, 
After  Williams,  Journal  of  Hydrology,  88:  pg  157) 
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Figure  3 


Relationships  between  Meander  Wavelength,  Channel  Width  and  Radius  of  Curvature 
over  a  wide  range  of  stream  size  { Rosgen,  D.,  Applied  River  Morphology  pg.  2-6;  After 
Leopold  1964) 


From  these  relationships  Williams^  derived  several  equations  of  best  fit  based  on  all  the 


data  available  to  him  at  the  time,  for  various  streams  and  rivers.  Additionally  an  equation 
(eq.7)  for  Thalweg  Depth  from  W.  B.  Leeder  is  also  given.  The  equations  of  interest  in 
the  estimation  of  depth  from  meander  geometry  measurements  are  as  follows: 

D=  0.12W**® 

eql. 

D  =  0.09  W*®**  K*** 

eq2. 

L»  =  240 

eq3. 

R«  =  42 

eq4. 

equations  3  and  4  can  be  rewritten  to  yield  dq>th  as  a  function  of  wavelength  and  radius  of  curvature 

E,  ^  p(lii(L,„/240)/1.52) 

oq5. 

jj  ^  ^(lii(Bc/42)/1.52) 

eq6. 

and  for  Thalweg  Depth 

W  =  «.*Dt 

eq7. 

which  can  aslso  be  rewritten  to  yield 

j,^^^(ln(ir/6.8)/1.54) 

Where 

D  =  (1ei)tlimeanliy(b'aii]k 
Dt  =  c1«i>th  thahwg 
W  =  width 

L]n  =  laeniulai' wnvelenglit 
=  radius  of  ctuvature 

K  =  sinuosity  (ititio  of  sti^nm  cliaiuiel  distanro  to  distance  down  slope) 


^  Williams,  G.P.,  1986.  River  meanders  and  channel  size.  J.  Hydrol.,  88: 147-164 
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Figure  8. 

Depth  as  a  function  of  Width  and  Sinuosity  (depth  curves  are  numbered  1  —  10  meters) 


In  order  to  help  the  reader  better  understand  the  terms  employed  in  the 
equations  the  Figures  9  and  10  are  provided. 


estimating 
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Flood  Depth 


Figure  9. 

Schematic  Illustration  of  the  different  “depth”  terms. 


Figure  10. 

Schematic  Illustration  of  Sinuosity 
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3.0  Method 


To  locally  verify  these  methods  of  depth  estimation  a  section  of  the  Wolf  River  located 
near  the  coast  of  southern  Mississippi  was  selected.  The  region  selected  for  testing  was 
assumed  to  be  sufficiently  inland  to  be  free  of  major  tidal  influence.  The  test  area  is 
shown  in  Figure  1 1 . 


Wolf  River  Test  Locations  -  River  Depth  Estimation 
Harrison  County  Mississippi 


552570  553570 


Figure  11. 

Wolf  River  test  locations. 
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USGS  Ortho-Photo  Quad  imagery  was  obtained  from  MARIS,  Mississippi  Automated 
Resource  Information  System,  in  the  Mississippi  Transverse  Mercator  projection.  All 
meander  measurements  were  taken  from  this  imagery  using  the  ESRI  ArcMap  GIS 
software.  For  each  of  the  five  test  locations,  river  width,  bend  radius,  meander 
wavelength  and  sinuosity  were  measured.  These  measurements  were  entered  into  the 
appropriate  equations  and  the  results  were  recorded.  To  facilitate  the  calculations  for 
repeated  usage,  a  Java  application  program  was  written.  This  program  is  portable  across 
a  wide  range  of  computer  operating  systems  and  is  given  in  the  appendix  of  this 
document.  In  addition  to  cdculating  a  depth  estimate  based  on  each  measured  quantity, 
this  program  also  computes  an  overall  estimate  based  on  an  average  of  each  of  the 
individual  calculations.  Since  the  measurement  of  meander  geometry  is  some  what 
subjective,  it  was  felt  that  an  average  of  the  four  methods  would  produce  more  consistent 
results.  These  averaged  values  were  used  in  the  analysis.  A  screen  capture  of  this 
application  is  shown  in  Figure  12. 


B  River  Depth  Estimator  vl2.1 1.2003 


Hie 


All  Units  Are  in  Meters 


Mean  Hydraulic  Depth 
Thalweg  Depth 

Mean  Hydaulic  Depth 


Meander  Wave  Length 


1 _ 
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Mean  Hydraulic  Depth 


Meander  Bend  Radius 


Figure  12.  Screen  Capture  of  the  Java  Application  User  Interface  used  for  calculating 
average  Mean  Hydraulic  Depth. 
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Field  measurements  were  also  made  at  each  of  the  test  locations.  Depth  cross-sections 
were  measured  using  a  lead-line  for  depth  and  “boat  lengths”  for  horizontal  displacement. 
Depth  accuracy  was  estimated  to  be  within  +/- 10  cm.  and  horizontal  distance  accuracy  to 
be  within  +/-  50  cm.  The  observed  cross-sections  were  plotted  and  mean  hydraulic 
depths  were  calculated  from  the  results.  An  area  for  each  cross  section  was  calculated  by 
trapezoidal  integration,  and  from  this  area  a  mean  depth  derived  that  when  multiplied  by 
the  measured  width  would  yield  an  equivalent  cross  sectional  area.  (i.e.  a  rectangle  of 
equal  area).  During  the  field  measurements  the  thalweg  depth'^  or  greatest  depth  value 
was  also  determined.  A  composite  of  the  cross-sections  observed  at  the  test  locations  is 
shown  in  Figure  13. 


ComposKe  X-Section 


“4—  Seriesi 
-c-“Series2 
SeriesS 
Series4 
“-X-"  S©ries5 


Figure  13. 

Composite  graph  of  the  measured  cross-sections  at  all  five  test  locations. 


thalweg  is  a  physical  geological  term  derived  from  two  German  words,  “Thai"  meaning  valley,  and 
“Weg”  meaning  path  or  trail,  thus  thalweg  mean  “  the  valley  path",  in  this  case  the  deepest  part  of  the 
river  channel. 
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4.0  Results 


Table  1  shows  the  results  for  the  five  selected  test  sites.  Width,  Wavelength,  Radius  of 
Curvature,  and  Sinuosity  as  measured  for  each  location  are  show  along  with  the 
calculated  and  field  measured  values  for  Mean  Hydraulic  Depth,  and  Thalweg  Depth. 
The  estimated  values  in  each  case  are  based  on  an  average  of  all  four  estimating 
equations,  i.e.  width,  wavelength,  radius  of  curvature  and  sinuosity  at  each  cross-section 
location.  The  expected  and  measured  results  were  in  close  agreement.  The  values 
associated  with  cross-section  #3  were  excluded  in  the  overall  error  estimate.  It  is  believe 
that  the  measured  Thalweg  Depth  for  this  location  was  located  behind  a  large  submerged 
log,  which  caused  excessive  bottom  scour  thus  making  the  measured  value  deeper  than 
normal.  The  columns  in  Table  1,  labeled  Measured  MHD  and  Measured  Thalweg,  and 
the  two  corresponding  Error  of  Estimation  columns  each  contain  two  sets  of  numbers. 
The  first  number  is  actual  water  depth  (float  a  boat)  and  the  second  number  is  an  adjusted 
depth  for  water  level  below  the  “Bank-full”  elevation.  It  is  the  second  number  in  this 
case  that  the  estimating  equations  are  designed  to  predict  rather  than  the  actual  water 
depth  observed.  However  the  actual  water  depth  is  also  given  as  the  first  number  of  the 
two. 


Sec 

# 

Width 

Wave 

Length 

Radius  of 
Curvature 

Sinuosity 

Estimated 

MHD 

Measured 

MHD 

Error  of 
Estimation 

Estimated 

Thalweg 

Measured 

Thalweg 

Error  of 
Estimation 

1 

40 

112 

2.1 

1.8 

3.1 

2.5  //2.8 

ItliV/lileH 

2 

139 

2.6 

1.9 

1.3// 1.6 

0.6//  0.3 

2.4 

2  //2.3 

3 

24 

PM 

52 

2.4 

1.4 

2.2//  2.5 

2.2 

3.6// 3.9 

”1.4  ** 

4 

31 

251 

55 

2.5 

1.5 

2.6 

2.5// 2.8 

0.1// -0.2 

5 

37 

\U8 

55 

2.9 

1.9 

1.6//  1.9 

0.3 //O 

2.9 

3  //3.3 

-0.1// -0.4 

All  units  are  in  Meters 

Bank  Full  Offset  of  0.3  meters 

**  values  not  used 

Mean  error 

0.25// 0.03 

.25// -0.05 

RMS 

0.35  //0.22 

Table  1,  Tabulated  Results 
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5.0  Conclusions 


The  method  employed  provided  a  reasonable  estimate  of  the  river  depths  in  the  areas 
selected.  Overall,  the  Root  Mean  Squared  Error  was  0.22  meters,  which  is  within  12 
percent  of  the  River  Mean  Hydraulic  Depth  in  the  test  area  and  an  RMS  of  0.27 
meters  for  the  thalweg,  which  is  within  10  percent  of  the  measured  depth. 

Additional  studies  are  needed  to  further  validate  this  method  for  a  wider  range  of  river 
types.  It  should  also  be  pointed  out  that  actual  water  depths  encountered  in  the  field  may 
vary  widely  from  these  estimates  as  a  result  of  climatic  trends  or  near  term  rain  fall 
amounts.  However,  the  estimating  techniques  are  designed  to  determine  the  “bank  full” 
depth  and  techniques  may  be  developed  to  adjust  these  estimated  values  based  on  a 
comparison  of  currently  observed  “river  width”  to  that  of  the  “bank  width”. 
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8.0  Appendix 

This  appendix  contains  the  Java  code  used  to  implement  the  depth  estimation  equations. 
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package  mil.navy.nrlssc.dmap.riverdepth; 

import  java.awt.*; 

import  java.awt.event.*; 

import  javax.swing.*; 

import  javax.swing.border.*; 

import  java.text.*; 

import  java.awt.print.*; 

import  j  avax.print.  * ; 

import  javax.print.attribute.*; 

import  javax.print.attribute.standard.*; 

import  javaio.*; 

import  java.awt.image.* ; 

import  javax.imageio.  * ; 

public  class  RiverDepthFrame  extends  JFrame  implements  ActionListener,  Printable  { 
JMenuBar  jmbMenu; 

JMenu  jmFile; 

JMenuItem  jmiPrint,  jmiSaveAs,  jmiExit; 

JLabeljlblTitle2, 

jlblBankWidthl,  jlblMeanHydraulicDepthl ,  jlblThalwegDepth, 
jlblMeanderW  aveLength,  jlblMeanHydraulicDepth2, 

jlblMeanderBendRadius,  jlblMeanHydraulicDepthS,  jlblBankWidth2, 
jlblSinuosity,  jlblMeanHydraulicDepth4,  jlblAverageMeanHydraulicDepth, 
JlblDepthRangeForAverageFlowConditions  1 ,  jlblComments; 

JTextField  jtfBankWidthl,  jtfMeanHydraulicDepthl,  jtfThalwegDepth, 
jtfMeanderWaveLength,  jtfMeanHydraulicDepth2, 

jtfMeanderBendRadius,  jtfMeanHydraulicDepthS,  jtfBankWidth2,  jtfSinuosity, 
jtfMeanHydraulicDepth4,  jtfAverageMeanHydraulicDepth; 

JTextArea  jtaComments; 

JButton  jbtnCalculatel ,  jbtnCalculate2,  jbtnCalculateS,  jbtnCalculate4,  jbtnAverage; 

DecimalFormat  oneFormatter; 

RiverDepthFrame  frameReference; 

PrintService  selectedPrinter; 

String  strVersion; 

JLabel  jlblPrintingVersion; 

String  strSaveReportAsImageFolder; 

public  RiverDepthFrame(){ 
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frameReference  =  this; 

setE)efaultCloseOperation(JFrame.EXIT_ON_CLOSE); 
setSize(620, 680); 

setLx)cation(  (int)(Tool]dt.getDefaultToolkit().getScreenSize().width/2.0  - 
getSize().width/2.0  +  0.5), 

(int)(Toolkit.getDefaultToolkit().getScreenSize().height/2.0  - 
getSize().height/2.0  +  0.5) ); 
setResizable(false) ; 
strVersion  =  "vl2. 12.2003"; 
setTitleC'River  Depth  Estimator "  +  strVersion); 

setIconImage(Utility.getImageFromJar("mil/navy/nrlssc/dmap/riverdepth/icons/Calculato 

r.gif')); 

getContentPaneO .  setLayout(null); 

UIManager.putC'Label.font",  new  Font( 
((Font)(UIManager.get("Label.font"))).getName(),  Font.PLAIN,  12)); 

UIManager.put("Button.font",  new  Font( 
((Font)(UIManager.get("Button.font"))).getName(),  Font.PLAIN,  12)); 

jmbMenu  =  new  JMenuBar(); 

jmFile  =  new  JMenu("File"); 
jmbMenu.add(jniFile); 

jmiPrint  =  new  JMenuItem(  "Print"); 

jiniPrint.addActionListener(this); 

jmFile.add(jniiPrint); 

jmiSaveAs  =  new  JMenuItem("Save  As..."); 

jmiSaveAs.addActionListener(this); 

jmFile.add(jmiSaveAs); 

jmiExit  =  new  JMenuItem("Exit"); 
jmiExit.addActionListener(this) ; 
jmFile.add(jniiExit); 

setJMenuBar(jmbMenu); 

jlblTitle2  =  new  JLabel("AIl  Units  Are  in  Meters",  SwingConstants.CENTER); 
jlblTitle2.setFont(  new  Font(  ((Font)(UIMariager.get("Label.font"))).getName(), 
Font.PLAIN,  20)); 

jlblTitle2.setSize(getSize().width  -  20,  25); 
jlblTitle2.setLocation(10, 10); 
jlblTitle2.setForeground(Color.red); 
getContentPane().add(jlblTitle2); 
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jtfBankWidthl  =  new  JTextField(); 
jtfBankWidthl.setSize(80,  25); 
jtfBankWidthl.setLocation(33, 70); 
getContentPane().add(j  tfBankWidth  1 ); 

jlblBankWidthl  =  new  JLabel("Bank  Width",  SwingConstants.CENTER); 
jlblBankWidthl.setSize(130, 25); 
jlblBankWidthl  .setLocation(10, 95); 
getContentPaneO.addOlblBankWidthl); 

jbtnCalculatel  =  new  JButton("Calculate"); 
jbtnCalculatel.setSize(80, 25); 
jbtnCalculate  1  .setLocation( 1 80, 70); 
jbtnCalculatel. setMargin(new  lnsets(0, 0,  0, 0)); 
jbtnCalculate  1  .addActionListener(this); 
getContentPane().add(jbtnCalculate  1 ); 

jtfMeanHydraulicDepthl  =  new  JTextField(); 
jtfMeanHydraulicDepth  1  .setSize(80, 25); 
jtfMeanHydraulicDepthl  .setLocation(330, 70); 
getContentPane().add(jtfMeanHydraulicDepth  1 ); 

jlblMeanHydraulicDepthl  =  new  JLabel("Mean  Hydraulic  Depth"); 
jlblMeanHydraulicDepthl  .setSize(  180, 25); 
jlblMeanHydraulicDepthl  .setLocation(420, 70); 
getContentPane().add(jlblMeanHydraulicDepthl); 

jtfThalwegDepth  =  new  JTextField(); 
jtfThalwegDepth.setSize(80, 25); 
jtfThalwegDepth.setLocation(330,  1 10); 
getContentPane().add(jtfrhalwegDepth); 

jlblThalwegDepth  =  new  JLabel("Thalweg  Depth"); 
jlblThalwegDepth.setSize(  1 80, 25); 
jlblThalwegDepth.setLocation(420, 1 10); 
getContentPane().add(jlblThalwegDepth); 

jtfMeanderWaveLength  =  new  JTextField(); 
jtfMeanderWaveLength.setSize(80, 25); 
jtfMeanderWaveLength.setLocation(33, 1 80); 
getContentPane().add(jtfMeanderWaveLength); 

jlblMeanderWaveLength  =  new  JLabel("Meander  Wave  Length", 
SwingConstants.CENTER); 

jlblMeanderW  aveLength.setSize( 1 30, 25) ; 
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jlblMeanderWaveLength.setLocation( 10, 205); 
getContentPane().add(jlblMeanderWaveLength); 

jbtnCalculate2  =  new  JButton("Calculate"); 
jbtnCalculate2.setSize(80, 25); 
jbtnCalculate2.setLocation(180,  180); 
jbtnCalculate2.setMargin(new  lnsets(0, 0, 0, 0)); 
jbtnCalculate2.addActionListener(this); 
getContentPaneO  .add(jbtnCalculate2) ; 

jtfMeanHydraulicDepth2  =  new  JTextField(); 
jtfMeanHydraulicDepth2.setSize(80, 25); 
jtfMeanHydraulicDepth2.setLocation(330, 180); 
getContentPane().add(jtfMeanHydraulicDepth2); 

jlblMeanHydraulicDepth2  =  new  JLabel("Mean  Hydaulic  Depth"); 
jlblMeanHydraulicDepth2.setSize( 1 80, 25); 
jlblMeanHydraulicDepth2.setLocation(420, 1 80); 
getContentPane().add(jlblMeanHydraulicDepth2); 

JtfMeanderBendRadius  =  new  JTextField(); 
jtfMeanderBendRadius.setSize(80, 25); 
jtfMeanderBendRadius.setLocation(33, 250); 
getContentPane().add(jtfMeanderBendRadius); 

jlblMeanderBendRadius  =  new  JLabel("Meander  Bend  Radius", 
SwingConstants.CENTER); 

jlblMeanderBendRadius.setSize(  130, 25); 
jlblMeanderBendRadius.setLocation(  10, 275); 
getContentPaneO  •  add(jlblMeanderBendRadius) ; 

jbtnCalculate3  =  new  JButton("Calculate"); 
jbtnCalculate3.setSize(80,  25); 
jbtnCalculate3.setLocation(180,  250); 
jbtnCalculate3.setMargin(new  lnsets(0, 0, 0, 0)); 
jbtnCalculate3  .add  ActionListener(this) ; 
getContentPaneO .  add(jbtnCalculate3) ; 

jtfMeanHydraulicDepth3  =  new  JTextFieldO; 
jtfMeanHydraulicDepth3.setSize(80,  25); 
jtfMeanHydraulicDepth3.setLocation(330,  250); 
getContentPaneO .  add(j  tfMeanHydraulicDepth3 ) ; 

jlblMeanHydraulicDepth3  =  new  JLabel("Mean  Hydraulic  Depth"); 
jlblMeanHydraulicDepth3.setSize(150, 25); 
jlblMeanHydraulicDepth3.setLocation(420, 250); 
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getContentPane().add(jlblMeanHydraulicDepth3); 

jtfBankWidth2  =  new  JTextFieldQ; 
jtfBankWidth2.setSize(80, 25); 
jtffiankWidth2.setLocation(33,  320); 
getContentPane().add(jtfBankWidth2); 

jlblBankWidth2  =  new  JLabel("Bank  Width",  SwingConstants.CENTER); 
jlblBankWidth2.setSize(130, 25); 
jlblBankWidth2.setLocation( 10, 345); 
getContentPane().add(jlblBankWidth2); 

jtfiSinuosity  =  new  JTextField(); 
jtfSinuosity.setSize(80, 25); 
jtfSinuosity.setLocation(33, 380); 
getContentPane().add(jtfSinuosity); 

jlblSinuosity  =  new  JLabel("Sinuosity",  SwingConstants.CENTER); 
jlblSinuosity.setSize(130, 25); 
jlblSinuosity.setLx)cation(  10, 405); 
getContentPane().add(jlblSinuosity); 

jbtnCalculate4  =  new  JButton("Calculate"); 
jbtnCalculate4.setSize(80, 25); 
jbtnCalculate4.setLocation(  1 80,  350); 
jbtnCalcuIate4.setMargin(new  lnsets(0, 0, 0, 0)); 
jbtnCalculate4.addActionListener(this); 
getContentPane().add(jbtnCalculate4); 

jtfMeanHydraulicDepth4  =  new  JTextFieldQ; 
jtfMeanHydraulicDepth4.setSize(80,  25); 
jtfMeanHydraulicDepth4.setLocation(330,  350); 
getContentPane().add(jtfMeanHydraulicDepth4); 

jlbIMeanHydraulicDepth4  =  new  JLabel("Mean  Hydraulic  Depth"); 
jlblMeanHydraulicDepth4.setSize( 1 80, 25); 
jlblMeanHydraulicDepth4.setLocation(420,  350); 
getContentPane().add5lblMeanHydraulicDepth4); 

jbtnAverage  =  new  JButton("  Average"); 
jbtnAverage.setSize(80, 25); 
jbtnAverage.setLocation(180, 440); 
jbtnAverage.setMargin(new  lnsets(0, 0, 0, 0)); 
jbtnAverage.addActionListener(this); 
getContentPane().add(jbtnAverage); 
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jtfAverageMeanHydraulicDepth  =  new  JTextField(); 
jtfAverageMeanHydraulicDepth.setSize(80, 25); 
jtfAverageMeanHydraulicDepth.setLcx:ation(330, 440); 
getContentPane().add(jtfAverageMeanHydraulicDepth); 

jlblAverageMeanHydraulicDepth  =  new  JLabel("Average  Mean  Hydraulic  Depth"); 
jlblAverageMeanHydraulicDepth.setSize(  1 80, 25); 
jlblAverageMeanHydraulicDepth.setLocation(420, 440); 
getContentPane().add(jlblAverageMeanHydraulicDepth); 

jlblComments  =  new  JLabel("Comments:"); 
jlblComments.setSize(80, 20); 
jlblCoimnents.setLocation(30, 500); 
getContentPaneO .  add(jlblComments) ; 

jtaComments  =  new  JTextArea(); 
jtaComnients.setSize(450, 105); 
jtaConinients.setLocation(110, 500); 
jtaCoinments.setBorder(new  LineBorder(Color.black)); 
getContentPane().add(jtaConmients); 

jlblPrintingVersion  =  new  JLabel("River  Depth  Estimator  "  +  strVersion, 

S  wingConstants.RIGHT) ; 

jlblPrintingVersion.setFont(new  Font( 
((Font)(UIManager.get("Button.font"))).getName(),  Font.ITALIC,  12)); 
jlblPrintingVersion.setSize(580,  20); 
jlblPrintingVersion.setLocation( 10,  660); 

oneFormatter  =  new  DecimalFormat("0.0"); 


//  connector  lines 

getContentPane().add(new  Line(jtfBankWidthl.getLocation().x  + 
jtfBankWidthl.getSizeQ.width,  jtfBankWidthl  .getLocationQ.y  + 
jtfBankWidth  1  .getSize().height/2, 

jbtnCalculatel.getLocation().x,  jtfBankWidth  l.getLocation().y  + 
JtfBankWidth  1  .getSize().height/2)); 

getContentPane().add(new  Line(jbtnCalculatel.getLocation().x  + 
jbtnCalculatel.getSize().width,  jbtnCalculatel.getLocation().y  + 
jbtnCalculate  1  .getSize().height/2, 

JtfMeanHydraulicDepth  1  !getLocation().x, 
jbtnCalculatel  .getLocation().y  +  JbtnCalculate  1  .getSize().height/2)); 

getContentPane().add(new  Line(JbtnCalculatel.getLocation().x  + 

JbtnCalculate  1  .getSize().width/2,  JbtnCalculate  1  .getLocation().y  + 

JbtnCalculate  1  .getSize().height, 
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jbtnCalculate  1  .getLocation().x  +  jbtnCalculatel  .getSize().  width/2, 
jtfThalwegDepth.getLocation().y  +  jtfThalwegDepth.getSize().height/2)); 

getContentPane().add(new  Line(jbtnCalculatel.getLocation().x  + 
jbtnCalculatel  .getSize().width/2,  jtfThalwegDepth.getLocation().y  + 
jtfThalwegDepth.getSize().height/2, 

JtfThalwegDepth.getLocationO.x, 

JtfThalwegDepth.getLocationO.y +jtfThalwegDepth.getSize().height/2)); 

getContentPane().add(new  Line(jtfMeanderWaveLength.getLocation().x  + 
JtfMeanderWaveLength.getSizeQ.width,  jtfMeanderWaveLength.getLocation().y  + 
jtfMeanderWaveLength.getSize().height/2, 

jbtnCalculate2.getLocation().x, 

JtfMeanderWaveLength.getLocationO.y  +  jtfMeanderWaveLength.getSize().height/2)); 

getContentPane().add(new  Line(jbtnCalculate2.getLocation().x  + 
jbtnCalculate2.getSize().width,  jbtnCalculate2.getLocation().y  + 
jbtnCalculate2.getSize().height/2, 

jtfMeanHydraulicDepth2.getLocation().x, 

jbtnCalculate2.getLocation().y+jbtnCalculate2.getSize().height/2)); 

getContentPane().add(new  Line(jtfMeanderBendRadius.getLocation().x  + 
jtfMeanderBendRadius.getSizeQ.width,  JtfMeanderBendRadius.getLocationO.y  + 
jtfMeanderBendRadius.getSize().height/2, 

jbtnCalculate3.getLocation().x, 

jtfMeanderBendRadius.getLx)cation().y  +  jtfMeanderBendRadius.getSize().height/2)); 

getContentPane().add(new  Line(jbtnCalculate3.getLocation().x  + 
jbtnCalculate3. gets  ize().  width,  jbtnCalculate3.getLocation().y  + 
jbtnCalculate3.getSize0.height/2, 

JtfMeanHydraulicDepth3.getLocation().x, 
jbtnCalculate3.getLocation().y  +jbtnCalculate3.getSize().height/2)); 

getContentPane().add(new  Line(jtfBankWidth2.getLocation().x  + 
jtfBankWidth2.getSize().width,  jtfBankWidth2.getLocation().y  + 
jtfBankWidth2.getSize().height/2, 

jbtnCalculate4.getLocation().x  -  30, 
jtfBankWidth2.getLocation().y  +  jtfBankWidth2.getSize().height/2)); 

getContentPane().add(new  Line(jtfSinuosity.getLocation().x  + 
JtfSinuosity.getSizeQ.width,  jtfSinuosity.getLocation().y  + 
jtfSinuosity.getSizeO.height/2, 

jbtnCalculate4.getLocation().x  -  30,  jtfSinuosity.getLocation().y  + 
jtfSinuosity.getSize().height/2)); 

getContentPane().add(new  Line(jbtnCalculate4.getLocation().x  -  30, 
jtfBankWidth2.getLocation().y  +  jtfBankWidth2.getSize().height/2, 

jbtnCalculate4.getLocation().x  -  30,  jtfSinuosity.getLocation().y  + 
jtfSinuosity.getSize().height/2)); 

getContentPane(Xadd(new  Line(jbtnCalculate4.getLocation().x  -  30, 
jbtnCalculate4.getLocation().y  +  jbtnCalculate4.getSize().height/2, 
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jbtnCalculate4.getLocation().x,  jbtnCalculate4.getLocation().y  + 
jbtnCalculate4.getSize().height/2)); 

getContentPane().add(new  Line(jbtnCalculate4.getLocation().x  + 
jbtnCalculate4.getSize().width,  jbtnCalculate4.getLocation().y  + 
jbtnCalculate4.getSize().height/2, 

jtfMeanHydraulicDepth4.getLocation().x, 
jbtnCalculate4.getLocation().y  +  jbtnCalculate4.getSize().height/2)); 

getContentPane().add(new  Line(jbtnAverage.getLocation().x  + 
jbtnAverage.getSizeQ.width,  jbtnAverage.getLocation().y  + 
jbtnAverage.getSize().height/2, 

jtfAverageMeanHydraulicDepth.getLocation().x, 
jbtnAverage.getLocation().y  +  jbtnAverage.getSize().height/2)); 

strSaveReportAsImageFolder  = 


} 

public  void  clearAverages(){ 
jtfAverageMeanHydraulicDepth.setTextC'"); 

} 

public  void  jmiExitClicked(){ 

System.exit(O); 

} 

public  void  jniiPrintClicked(){ 

(new  PrintDialog(this)).setVisible(true); 

} 

class  ImageWriterThread  extends  Thread{ 

File  outputFile; 

String  strlmageType; 

public  ImageWriterThread(File  outputFile,  String  strlmageType)  { 
this.outputFile  =  outputFile; 
this.strLtnageType  =  strlmageType; 

} 

public  void  run(){ 
try{ 

setCursor(new  Cursor(Cursor.WAlT_CURSOR)); 

Bufferedimage  panellmage  =  new 

BufferedImage(getContentPane().getSize().width,  getContentPane().getSize().height, 
BufferedImage.TYPE_INT_RGB); 

Graphics2D  panelImageG2D  =  panelImage.createGraphics(); 
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jlblPrintingVersion.setLocation(jlblPrintingVersion.getLocation().x, 
getContentPane().getSize().height  -  20); 

getContentPaneO  .add(j  IblPrinting  Version) ; 

getContentPane().paint(panelImageG2D); 

getContentPane().remove(jlblPrintingVersion); 

ImageIO.write(panelImage,  strlmageType,  outputFile); 
JOptionPane.showMessageDialog(frameReference,  "Report  saved  as  "  + 
outputFile,  "Save  Successful",  JOptionPane.INFORMATION_MESSAGE); 

} 

catch(Exception  e){ 

JOptionPane.showMessageDialog(frameReference,  "Error  Avriting  file: "  + 
e.getMessageO, 

"Save  Failed",  JOptionPane.ERROR_MESSAGE); 

} 

setCursor(new  Cursor(Cursor.DEFAULT_CURSOR)); 

} 

} 

public  void  jmiSaveAsClicked(){ 
frameReference.requestFocusO; 

JFileChooser  saveReportAsFileChooser  =  new  JFileChooser(); 
saveReportAsFileChooser.setCurrentDirectory(new 
File(strSaveReportAsImageFolder)); 

saveReportAsFileChooser.setDialogTitleC'Save  As..."); 
saveReportAsFileChooser.setAcceptAllFileFilterUsed(false); 

saveReportAsFileChooser.addChoosableFileFilter(new  GenericFileFilter(".png", 
"PNG  Files")); 

saveReportAsFileChooser.addChoosableFileFilter(newGenericFileFilter(".jpg", 

"JPG  Files")); 

saveReportAsFileChooser.setFileFilter((saveReportAsFileChooser.getChoosableFileFilte 

rs())[0]); 

int  approveSaveValue; 

int  approveOverwriteValue  =  JOptionPane.NO_OPTION; 
while(!  (approveOverwriteValue  =  JOptionPane.YES_OPTION)){ 
approveSaveValue  =  saveReportAsFileChooser.showSaveDialog(tliis); 
if  (approveSaveValue  !=  JFileChooser. APPROVE_OPTION){ 
return; 

} 

if  ( saveReportAsFileChooser.getSelectedFile().exists()  II 

(new  File(saveReportAsFileChooser.getSelectedFile().toString()  + 
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((GenericFileFilter)(saveReportAsFileChooser.getFileFilter())).getExtension())).exists() 
){  //  show  "replace?"  dialog 

approveOverwrite Value  =  JOptionPane.showConfimiDialog(this,  "File  already 
exists,  overwrite?",  "Confirm  Overwrite", 

JOptionPane.YES_NO_OPTION, 

JOptionPane.WARNING.MESSAGE); 

} 

else{ 

approveOverwriteValue  =  JOptionPane.YES_OPTION; 

} 

} 

strSaveReportAsImageFolder  = 
saveReportAsFileChooser.getCurrentDirectoryO.toStringO; 

String  strSelectedFileFilterExtension  = 

((GenericFileFilter)(saveReportAsFileChooser.getFileFilter())).getExtension(); 

String  strSelectedFile  =  saveReportAsFileChooser.getSelectedFile().getPath(); 
if(! 

(strSelectedFile.toLowerCase().endsWith(strSelectedFileFilterExtension.toLowerCase())) 

){ 

StrSelectedFile  =  strSelectedFile  +  strSelectedFileFilterExtension; 

} 

ImageWriterThread  imageSaver  =  new  ImageWriterThread(new  File(strSelectedFile), 
strSelectedFileFilterExtension.substring(  1 )); 
imageSaver.  start() ; 

} 

public  void  jbtnCalculate  1  Clicked()  { 
clearAveragesO; 
try{ 

jtfMeanHydraulicDepth  1  .setText(oneFonnatter.format(0. 12  * 
Math.pow(Double.parseDouble(jtfBankWidth  1  .getTextQ),  0.69))); 

jtfrhalwegDepth.setText(oneFormatter.format(0.28  * 
Math.pow(Double.parseDouble(jtfBankWidth  1  .getTextO),  0.65))); 

} 

catch(Exception  e){ 

JOptionPane.showMessageDialog(this,  "Please  input  a  valid  number.",  "Invalid 
Number",  JOptionPane.ERROR_MESSAGE); 

} 

} 

public  void  jbtnCalculate2Clicked()  { 
clearAveragesO; 
try{ 
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jtfMeanHydraulicDepth2.setText(oneFormatter.format(Math.exp(Math.log(Double.parse 
Double(jtfMeanderWaveLength.getText())  /  240.0)  /  1.52))); 

} 

catch(Exception  e){ 

JOptionPane.showMessageDialog(this,  "Please  input  a  valid  number.",  "Invalid 
Number",  JOptionPane.ERROR_MESSAGE); 

} 


} 

public  void  jbtnCalculate3Clicked(){ 
clearAveragesO; 
try{ 

jtfMeanHydraulicDepth3.setText(oneFormatter.format(Math.exp(Math.log(Double.parse 
Double(jtfMeanderBendRadius.getText())  /  42.0)  /  1.52))); 

} 

catch(Exception  e){ 

JOptionPane.showMessageDialog(this,  "Please  input  a  valid  number.",  "Invalid 
Number",  JOptionPane.ERROR_MESSAGE); 

} 

} 

public  void  jbtnCalculate4Clicked()  { 
clearAveragesO; 
try{ 

jtfMeanHydraulicDepth4.setText(oneFormatter.format(0.09  * 
Math.pow(Double.parseDouble(jtffiankWidth2.getText()),  0.59)  * 

Math.pow(Double.parseDouble(jtfSinuosity.getText()),  1 .46))); 

} 

catch(Exception  e){ 

JOptionPane.showMessageDialog(this,  "Please  input  valid  numbers.",  "Invalid 
Number",  JOptionPane.ERROR_MESSAGE); 

} 

} 

public  void  jbtnAverageClicked(){ 
try{ 

jtfAverageMeanHydraulicDepth.setText(oneFormatter.format( 
(Double.parseDouble(jtfMeanHydraulicDepth  1  .getText())  + 

Double.parseDouble(jtfMeanHydraulicDepth2.getText())  + 

Double.parseDouble(jtfMeanHydraulicDepth3.getText())  + 
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Double.parseDouble(jtfMeanHy(iraulicDepth4.getText()))  /  4,0 )); 

} 

catch(Exception  e){ 

JOptionPane.showMessageDialog(this,  "All  4  values  not  yet  calculated.",  "Missing 
Calculation(s)",  JOptionPane.ERROR_MESSAGE); 
return; 

} 

} 

public  void  actionPerformed(ActionEvent  ae){ 

if  (ae.getSourceO  ==  jiniPrint){ 
jmiPrintClickedO; 

} 

else  if  (ae.getSourceO  ==  jniiSaveAs){ 
jmiS  ave  AsClickedO ; 

} 

else  if  (ae.getSourceO  ==  jmiExit){ 
jmiExitClickedO; 

} 

else  if  (ae.getSourceO  ==  jbtnCalculatel){ 
jbtnCalculate  1  Clicked(); 

} 

else  if  (ae.getSourceO  ==  jbtnCalculate2){ 
jbtnCalculate2Clicked(); 

} 

else  if  (ae.getSourceO  ==  jbtnCalculate3){ 
jbtnCalculate3Clicked(); 

} 

else  if  (ae.getSourceO  ==  jbtnCalculate4){ 
jbtnCalculate4Clicked(); 

} 

else  if  (ae.getSourceO  ==  jbtnAverage){ 
jbtn  AverageClickedO ; 

} 

} 

class  Line  extends  JPanel{ 

public  Line(int  startX,  int  start Y,  int  endX,  int  endY){ 
setLocation(startX,  start Y); 
setSize(endX  -  startX  +  1,  endY  -  startY  +1); 
setBackground(Color.black); 

} 

} 

public  static  void  main(String[]  args){ 
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RiverDepthFrame  riverDepthScreen  =  new  RiverDepthFrame(); 
riverDepthScreen.setVisible(true); 

} 


public  int  print(Graphics  g,  PageFonnat  pf,  int  pageindex) 
throws  PrinterException  { 
if  (pageindex  >=  1)  { 
return  Printable.NO_SUCH_PAGE; 

} 

double  pageHeight  =  pf.getlmageableHeightO; 
double  pageWidth  =  pf.getImageableWidth(); 


double  widthScaleFactor  =  pageWidth/(double)(getContentPane().getSize().width); 
double  heightScaleFactor  =  pageHeight/(double)(getContentPane().getSize().height); 

((Graphics2D)g).translate(pf.getImageableX(),  pf.getImageableYO); 
((Graphics2D)g).scale(Math.min(widthScaleFactor,  heightScaleFactor), 
Math.min(widthScaleFactor,  heightScaleFactor)); 

RepaintManager  currentManager  = 
RepaintManager.currentManager(getContentPane()); 
currentManager.setDoubleBufferingEnabled(false); 

getContentPane().paint(g); 

currentManager.  setDoubleBufferingEnabled(true) ; 

return  Printable.PAGE_EXISTS; 

} 

class  PrintThread  extends  Thread  { 
public  void  run(){ 
try{ 

setCursor(new  Cursor(Cursor.WArr_CURSOR)); 
frameReference.requestFocusO; 

PrinterJob  printJob  =  PrinterJob.getPrinterJob(); 

PageFormat  pf  =  new  PageFormat(); 

printJob .  setPrintable(frameReference,  pf) ; 
printJob.setPrintService(selectedPrinter); 

PrintRequestAttributeSet  printSettings  =  new  HashPrintRequestAttributeSet();  // 
set  up  printing  attributes 

printSettings.add(new  Copies(  1 )); 

jlblPrintingVersion.setLocation(jlblPrintingVersion.getLocation().x, 
getContentPane().getSize().height  -  20); 
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getContentPaneO .  add(j  IblPrinting  Version) ; 

printJob.print(printSettings); 

getContentPane().remove(jlblPrintingVersion); 

} 

catch(Exception  e){ 
e.printStackTraceO; 

} 

setCursor(new  Cursor(Cursor.DEFAULT_CURSOR)); 

} 

} 


class  PrintDialog  extends  JDialog  implements  ActionListener{ 

JLabel  jlblPrinter; 

JComboBox  jcbPrinters; 

JButton  jbtnOK,  jbtnCancel; 

PrintService[]  availablePrinters; 
public  PrintDialog(JFrame  parent)  { 
super(parent,  "Select  Printer",  true); 
setSize(400, 130); 

setLocation(  (int)(Toolkit.getDefaultToolkit().getScreenSize().width/2.0  - 
getSize().width/2.0  +  0.5), 

(int)(Toolkit.getDefaultToolkit().getScreenSize().height/2.0  - 
getSize().height/2.0  +  0.5) ); 

getContentPane().setLayout(null); 

jlblPrinter  =  new  JLabel("Printer:"); 
jlblPrinter.setSize(80, 20); 
jlblPrinter.setLocation(10, 10); 
getContentPane().add(jlblPrinter); 

jcbPrinters  =  new  JComboBox(); 
jcbPrinters.setSize(380,  20); 
jcbPrinters.setLocation(10,  30); 
getContentPane().add(jcbPrinters); 

jbtnOK  =  new  JButton("OK"); 
jbtnOK.setSize(50, 25); 
jbtnOK.setLocation(140, 70); 
jbtnOK.setMargin(new  lnsets(0, 0, 0, 0)); 
jbtnOK.addActionListener(this); 
getContentPane().add(jbtnOK); 

jbtnCancel  =  new  JButton("Cancel"); 
jbtnCancel.setSize(50, 25); 
jbtnCancel.setLocation(200, 70); 
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jbtnCancel.setMargin(new  lnsets(0, 0, 0, 0)); 
jbtnCancel.addActionListener(this); 
getContentPane().add(jbtnCaiicel) ; 

frameReference.setCursor(new  Cursor(Cursor.  W  A]T_CURS  OR)) ; 
availablePrinters  =  PrinterJob.lookupPrintServices(); 
frameReference.setCursor(new  Cursor(Cursor.DEFAULT_CURSOR)); 
Object  selectedPrinter  =  jcbPrinters.getSelectedItem(); 
jcbPrinters.remove  AllltemsO ; 
for  (int  i  =  0;  i  <  availablePrinters.length;  i++){ 
jcbPriiiters.addItem(availablePrinters[i].getName()); 
if  ( (selectedPrinter  !=  null)  &&  ( 

((String)(selectedPrinter)).equals(availablePrinters[i].getName())) ){ 
jcbPrinters.setSelectedlndex(i); 

} 

} 

} 

public  void  jbtnOKClicked(){ 

selectedPrinter  =  availablePrinters[jcbPrinters.getSelectedIndex()] ; 

(new  PrintThread()).start(); 
setVisible(false); 

} 

public  void  jbtnCancelClicked(){ 
setVisible(false); 

} 

public  void  actionPerformed(ActionEvent  ae){ 
if  (ae.getSourceO  ==  jbtnOK){ 
jbtnOKClickedO; 

} 

else  if  (ae.getSourceO  ==  jbtnCancel){ 
jbtnCancelClickedO; 

} 

} 

} 


} 


package  mil.navy.nrlssc.dmap.riverdepth; 

import  javaJo.*; 

public  class  StreamMethods{ 


public  static  byte[]  fullyReadStream(InputStream  streamln)  throws  Exception  { 
byte[]  bytesInStream  =  new  byte[524288];  //  initially  512  k 
byte[]  newBytesInStream; 
byte[]  buffer  =  new  byte[65536]; 
int  bytesReadThisPass  =  0; 
int  totalBytesRead  =  0; 

while  ( (bytesReadThisPass  =  streainln.read(buffer,  0,  buffer.length))  !=  -1){ 
while  ( (bytesReadThisPass  +  totalBytesRead)  >  bytesInStreamdength ){ 
newBytesInStream  =  new  byte[bytesInStream.length  +  524288]; 
System.arraycopy(bytesInStream,  0,  newBytesInStream,  0,  totalBytesRead); 
bytesInStream  =  newBytesInStream; 

} 

System.arraycopy(buffer,  0,  bytesInStream,  totalBytesRead,  bytesReadThisPass); 
totalB54esRead  =  totalBytesRead  +  bjdesReadThisPass; 

} 

newBytesInStream  =  new  byte[totalB5desRead]; 

System.arraycopy(bytesInStream,  0,  newBytesInStream,  0,  totalBytesRead); 
return  newB5desInStream; 

} 

} 
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package  mil.navy.nrlssc.dmap.riverdepth; 


import  java.io.File; 

import  javax.swing.*; 

import  j  avax.swing.filechooser.  * ; 

public  class  GenericFileFilter  extends  FileFilter{ 

String  strExtension,  strDescription; 

public  GenericFileFilter(String  strExtension,  String  strDescription)  { 
this.  strExtension  =  strExtension; 
this.strDescription  =  strDescription; 

} 

public  boolean  accept(File  fileToCheck){ 
if  (fileToCheck.isDirectory())  { 
return  true; 

} 

if  (fileToCheck.getName().toLowerCase().endsWith(strExtension.toLowerCase())){ 
return  true; 

} 

return  false; 

} 

public  String  getExtension(){ 
return  strExtension; 

} 

public  String  getDescriptionQf 
return  strDescription  +  "  (*"  +  strExtension  +  ")"; 

} 


} 
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package  mil.navy.nrlssc.dmap.riverdepth; 


import  java.awt.*; 
import  java.io.*; 
import  javax.swing.*; 

public  class  Utility  { 

public  static  Image  getImageFromJar(String  strImagePathInJar){ 

Image  loadedimage  =  null; 
try{ 

InputStream  imageStream  = 

Thread.currentThread().getContextClassLoader().getResourceAsStream(strImagePathInJ 

ar); 

loadedimage  = 

Toolkit.getDefaultToolkit().createImage(StreamMethods.fullyReadStream(imageStream) 

); 

imageStream.closeO; 

MediaTracker  mt  =  new  MediaTracker(new  JPanelO); 

mt.addImage(loadedImage,  0); 

mt.waitForAll(); 

} 

catch(Exception  e){ 
e.printStackTraceO ; 

} 

return  loadedimage; 

} 

} 
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